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Abstract 

As a generalization of a previous work [Phys. Rev. D. 59, 105014 (1999)], 
we compute analytically a class of three-loop vacuum diagrams with two ar- 
bitrarily different mass scales. We use a decomposition algorithm in which 
the integrand of the final integral for the third momentum vector, say, k, 
becomes independent of the angles of A;- vector in spherical polar coordinates. 
This algorithm proves to be very efficient in obtaining exclusively all e-pole 
terms of the given diagram. 
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Vacuum diagrams appear in the evaluation of the effective potential as well as in low- 
energy limits of certain physical amplitudes or Taylor coefficients of multipoint Green func- 
tions. After the algebra of Dirac gamma matrices and Lorentz indices is performed, any 
vacuum diagram can be reduced to some linear combinations of scalar vacuum integrals 
By the Chetyrkin-Misiak-Miinz algorithm 0, the computation of even non- vacuum 
diagrams can be reduced to the problem of the scalar vacuum diagrams. 

The authors have computed, using the Avdeev algorithm of recurrence relations [|T|, 
the single- (mass-) scale three-loop vacuum diagrams appearing in the three-loop effective 
potential calculation of the massive 0^ theory to the order of the finite term in the e expansion 
0. In the minimal subtraction scheme, knowledge of the finite terms as well as the pole terms 
is needed in the renormalization of the effective potential P-^]. However, in nonminimal 
subtraction schemes, the so-called 'renormalization conditions' are imposed on the effective 
potential • In these schemes, it is sufficient to know only the pole terms of the diagrams. 

It is very difficult to apply the Avdeev algorithm, which was written for single-scale 
three-loop vacuum diagrams, to the computation of two-scale three-loop diagrams. For the 
effective potentials of the 0{N) cj)'^ theory, it is necessary to compute two-scale three- loop 
vacuum diagrams [|1^]. Note that for the massive 0{N) 0^ theory, the mass-squared ratio 
of the two scales is arbitrary whereas for the case of the massless 0{N) cj)'^ theory, the ratio 
is 1/3. 

In a previous paper p|, we computed a class of two-scale three- loop vacuum diagrams 
with mass-squared ratio of 1/3 by applying the Chetyrkin-Misiak-Miinz algorithm, which 
was written also for a single-scale three-loop vacuum diagrams, to our two-scale problems. 
The purpose of this work is to evaluate analytically three-loop vacuum diagrams with two 
arbitrarily different mass scales by use of the Chetyrkin-Misiak-Miinz algorithm. In this 
algorithm, the integrand of the final integral for the third momentum vector, say, k, becomes 
independent of the angles of fc-vector in spherical polar coordinates. This algorithm proves 
to be very efficient in obtaining exclusively all e-pole terms of the given diagram. 

The integrals we consider are the following seven nontrivial vacuum integrals with two 
mass scales (m^ and me): 
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where Q{k'^;x) {x = m\^m\) and Ti^k"^ ; m\) are propagator-type one-loop integrals: 
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Throughout the paper, the momenta appearing in the formulas are all (Wick-rotated) Eu- 
clidean ones. The abbreviated integration measure is defined as 
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where D = 4 — e is the space-time dimension in the framework of dimensional regularization 
and /i is an arbitrary constant with mass dimension. Notice that the routing of (internal) 
momenta for J{h) is possible in either way given in Eq. (|l]). 

The functions Q{k^;x) and Ti.{k'^]m\,m\) in Eq. (^ have simple poles in e-} 
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With these splittings, we decompose all of the integrals in Eq. (|T|) as follows: 
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^Although the exact forms of G{k'^;x) and H (k'^ ; m\, ml) are available |11|, it is a very difficult 
task to carry out the /^-integrations in Eq. (0) below exactly. 
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Here, 5*1 to ^4 and Wi to Ws are lower-loop integrals [Q, whose definitions and values are 
relegated to the Appendix. In an abuse of terminology, we will call each first term of the 
right-hand sides of Eq. (^) a 'genuine' three-loop integral for the corresponding three-loop 
vacuum diagram. A remarkable feature of the above decompositions is that the integrands 
of genuine three-loop integrals are, in spheical polar coordinates, independent of the angle 
variables of the (4 — e) -dimensional fc-vector. 

A finite volume integration of each k integral in Eq. cannot give any e-pole term. 
Such a pole term can only arise from integration over large fc^ in Eq. (^. Therefore, knowing 
the behaviors of x) and H{k^; m\, m%) at large is enough to find the UV-divergent 
parts of all of the integrals in Eq. (p. The two-scale propagator-type one-loop inegral, 
m\, ml), can be expanded as an infinite sum of single-scale propagator-type one-loop 
integrals: 
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The propagator-type one-loop integral with an arbitrary integer power of the second prop- 
agator was calculated in a closed form |1TI| : 
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In order to see the asymptotic behaviors of G{k'^;x) and H{k'^;m\,m%), we must expand 
the hypergeometric function 3F2 in Eq. (j^) in an asymptoic form. We simply quote the 
result 0: 
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Note that a''''"(l,0) does contain a simple pole which cancels the pole of 6(1,0) in the 
expression F^'j^{k'^). Only a few lowest terms in the expansion of Eq. (|^ affect the pole 
parts of the considered three-loop integrals. 

The case of a = 1 in Eq. gives Q{k^, m\). Thus, the asymptotic behavior of G{k^; m\) 
is given by 
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From this G{k'^\ ^"a), the asymptotic behavior of G{k'^\ m^) can be written as 
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From Eqs. (^ and (^, we obtain the asymptotic behavior of if(A;^; m^, m^) in a double 
series form: 
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We are now in a position to carry out the k integrations in Eq. 
the following integral: 
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By doing first a trivial angle integration in (4— e) -dimensional spherical polar coordinates and 
then changing the remaining integration variable k, the magnitude of the (4 — e) -dimensional 
momentum, into m^/k"^, one can arrive at 
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where the symbol (~") is the expansion coefficient of (1 -|- = Yl'^o{~P)x''. As was 
mentioned above, the finite-volume k integration does not give rise to any e-pole term. 
Thus, within the desired accuracy of e-pole terms, we discard the finite-volume integration, 
the second integral J^yk'^ d{m\/A'^) ■ ■ ■. From the result of the trivial integration for a single 
variable m\/ k"^ — trivial in the sense that the integrand is a series sum of power functions of 
m\/k^ — it is not difficult to sort out pole terms: 
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Prom this equation, three integrals corresponding to n = 0, 1, 2 are evaluated as follows: 



(47r)H47r//2 



m 



B 



1 

+ - 
e 



V47r//2 

64 



-c/2 



/ G{k';ml)G{k'; 

Jk 



nil 



26 / _ m| 
3eH 



m. 



(f-Hi-(-T-^)(-(i)><-i)'<3)} 



80 



10 
"3 

27r2 



- 327 + 127" + ^)-f+(-7^+7-7'-7r 



11 



TT 



6 



6 



2 

£ 

1 



.(4-47.(-16.127):|.(l-27)(;|)>n(;|) 



i 



-(1-3) ^A'^]\+o{<P) 



m. 



mz 



-c/2 



(47r)^ V47r/x2/ \47r/x 

1 rl6 



(m^y/^ r G{e;ml)G{k^;ml) 



k"^ + m\ 



8 16 m| 
3 3 m\ 



1 flS , , .rni, ( ,mi^\ ^ /^rM 

+ _{__4,+ (4-4,)-|-(2 + 4-|)l„(-|)} 



TT 



+ -<^-2 + 7" + — + 



6 



h 7 

3 ' 



6 / 



i 

■7 



mi 



mi 



m 



47 + 27^ In 



(4) 



6 



(47r)4 V47r/i 



+ 



+ In 



1 



3e3 



2 



3 

+ 0(e°) 



In the above, the overall multiplying factors are defined as 
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and 7 is the usual Euler constant, 7 = 0.5772156649 ■ ■ ■. 

With Eq. (|l3), in a manner similar to that leading to Eq. (|Tl|) , we arrive at 
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from which the three integrals for n = 0,1,2 are calculated as follows: 
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Finally, using the asymptotic form of G{k'^; rn^), Eq. (^, we obtain 
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After substitutions of the appropriate values for and hj., we end up with 
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In summary, we tabulate a class of three-loop vacuum integrals given in Eq. (|Tp by 
substituting the results for the 'genuine' three-loop integrals scattered in Eqs. (0), (|14|), 
and (|T^) and for the 'factorized' lower-loop integrals in Eq. ( [AID i'^^o Eq. (^. Since, in 
Eq. (^), one- loop integrals {Si to 5*4) are multiplied by 1/e^ and two- loop integrals (Wi to 
M/g) by 1/e, it is sufficient to know their values to the order of and e°, respectively. The 
tabulation reads as follows: 
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APPENDIX: LOWER-LOOP INTEGRALS 



We give the definitions of lower-loop integrals Si to and Wi to Ws appearing in Eq. 
and their values: 
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